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Abstract 

We show that quark models of current matrix-elements based on the Bakamjian- 
Thomas construction of relativistic states with a fixed number of particles, plus 
the additivity assumption, are covariant in the heavy-quark limit and satisfy 
the full set of heavy-quark symmetry relations discovered by Isgur and Wise. 
We find the lower bound of p 2 in such models to be 3/4 for ground state mesons, 
independently of any parameter. Another welcome property of these models 
is that in the infinite momentum limit the wave functions vanish outside the 
domain < x < 1. 
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1 Introduction 



The necessity of a relativistic treatment of hadron center-of-mass motion in quark 
models is manifest in the calculation of form factors at high three-dimensional mo- 
mentum transfer q and effort in this direction have been made since a long time. 
On the other hand, it has been realized rather recently that in QCD hadronic form 
factors must satisfy a set of remarkable relations in the limit where the mass of the 
active quark in the hadron is made heavy [jlj, the so-called heavy quark symmetries. 

In the past few years ago 0, we have noticed that models based on a very 
simple treatment of hadron motion with Lorentz boost of spins and Lorentz contrac- 
tion of spatial wave functions - which we have formulated in the 70's || - indeed 
present these symmetry properties in the heavy quark limit 0, ||. This is not at 
all trivial because, as can be observed in the literature, most current models do not 
satisfy these properties [Q (or enforce them by hand). However our own models Q 
M have serious drawbacks. They are not covariant and they use, in addition to the 
basic assumptions of quark models, a series of approximations which are not settled 
in a well-defined framework. In addition, they do not show the expected behavior 
at large q 2 , which in turn is related to the fact that the null-plane limit of the wave 
functions do not vanish outside the domain < x < 1. 

We will use here such a well-defined framework which in addition turns out to 
have the outstanding merit of being covariant in the heavy quark limit and maintain 
the scaling properties found in the naive models. The purpose of this letter is to show 
these properties of this new approach and to derive a bound on the widely discussed 
parameter p 2 , the slope of the Isgur-Wise function. In addition, the class of models 
obtained in that way solves the above-mentioned problem at large q 2 , since, in the 
infinite momentum limit, the wave functions vanish outside the domain < x < 1, 
as will be proved elsewhere. 

Let us say, before proceeding, a few general words on the method. Although 
this does not seem to be expected from the knowledge of field theory, it is a very 
old finding || , , |§ that one can make an important step towards a fully relativistic 
theory with a fixed number of interacting constituents, with wave functions imple- 
menting a representation of the Poincare group through the construction of the full 
set of generators, and with a rest-frame Hamiltonian (or mass operator) containing 
rather standard (non-relativistic looking) potentials depending on relative coordi- 
nates, two-dimensional Pauli spins, . . . 

In this framework, the problem of knowing the relativistic wave function in 
motion in terms of the wave function at rest is solved exactly and in a rather simple 
manner, through a change of variables which is known once for all and does not 
depend on the interaction. 

The whole solution relies on a complete separation between two types of vari- 
ables, related to individual particle variables by explicit expressions which are simple 
at least in momentum space. On the one hand, we have global variables which 
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describe the whole system in analogy with the one-particle state. On the other 
hand, we have internal variables which somewhat generalise the relative variables of 
non-relativistic systems, a major property being that the two types of variables are 
commuting as operators. 

The important point is then that one can construct the Poincare generators as 
the ones of a free particle described by the global variables (total momentum P and 
relativistic center-of-mass position R) and a mass which can be chosen arbitrarily 
provided it depends only on the internal variables. 

This framework would be ideally suited to formulate relativistic quark models 
except for one serious drawback. There is no known covariant current operator, 
and in particular the usual one-quark free current, which corresponds to the basic 
additivity assumption of quarks models, is not covariant in general when sandwiched 
between hadron states. This failure entails that we have not a satisfactory relativistic 
model for transitions, which was however precisely the initial motivation for appealing 
to this treatment. We have then to return in general once more to the old discussions 
on a possible approximate covariance under certain particular conditions, or on the 
choice of a best reference frame. 

One remarkable exception to this failure is the case of the heavy quark limit. 
Let us recall that this limit consists in considering systems containing one quark very 
heavy with respect to all the others, and moreover to consider transitions between 
such "heavy-light" hadrons, where it is the heavy quark which endows electroweak 
interaction through some external current. In this limit, we find that the above 
framework with the one-quark current gives a covariant model for transitions. It 
is then on this limit that we shall concentrate after a presentation of the general 
formalism, leaving the study of other situations for future discussions. We recall that 
the interest of this limit is not purely theoretical. It is believed to be roughly realized 
in the B -> D^lv semi-leptonic decays, where one is measuring in particular with 
increasing accuracy the slope p 2 at the no-recoil point, and where the knowledge of 
form factors is expected to improve much in the future. It is interesting that our 
lower bound on p 2 is rather close to the value estimated by CLEO II experiment. 

On the theoretical side, we must emphasize that the model presented in this 
letter is not just one more quark model. It has the interest of embedding many 
recent attempts. Indeed, we find that certain recent models for mesons or baryons 
are explicitly based on the B-T formalism plus the free quark current f5||, |jl3[| . On the 
other hand, there are many approaches directly formulated on the null-plane, among 
which several \ TU \, ]T2|] , can be shown to be the P = oo limit of the present 



approach. They will have the same heavy quark limit obtained in this letter, since 
we show this limit to be covariant. Also, according to our findings, it seems that 
the intuitive approach of Close and Wambach directly formulated for the heavy 
quark limit, leads to the same final expressions. One consequence is that the lower 
bound on p 2 we have found is of general interest in that it will apply to a large class 
of proposed models. 
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Of course, this does not apply to all models encountered in the literature, and 
one may find seemingly similar approaches which differ by certain details which, in 
the end, reveal crucial \TA\. We have also definitely different approaches leading to 



covariant and scaling models, like the one of Kaidalov |T6 

As to our previous model 0, it comes out that it could be defined, at the start, 
as an approximation to the present one, which would consist mainly in neglecting 
internal momenta with respect to the quark masses. However, although this could 
seem quite natural in a quark model, it is found to lead finally to drastic differences 
with the present one, especially at large q 2 , on which we will comment elsewhere. 



2 Relativistic quark models for currents in the 
Bakamjian-Thomas formalism. 

In the present section, we will not yet take the infinite mass limit. The n-particle 
Hilbert space, which is naturally the tensor product of the n individual one-particle 
Hilbert spaces, is made of functions Sl ,...s n (pu ■ ■ -Pn) °f the so-called one-particle 
variables, spins S{ and momenta pi. Assuming that particle 1 is the active particle, 
the additivity hypothesis means that the current density operator in the n-particle 
Hilbert space will be the tensor product of the current density operator on particle 1 
by the identity operators on particles 2, . . . n. This writes : 

(V\0\V) = / dfl . A (f[ - jj* ) T V (l) 

S ' ' ' J (2nf(2n) 3{ lU^) 3} fa s fa Sn 1 ) 

K' v s 2 ,...s n (Pl,P2, ■ ■ ■Pn)*0(pi,Pl) s > l , si y suS2 ,... Sn (p 1 ,p^ ...p n ) 

where from now on the primes will denote the final states, and 0(p',p) s i jS is the 
matrix element between one-particle states: 

0(p\p) s ,, s = (p\s'\0\p } s) (2) 

The one-particle states that we use are defined, including their normalisation, by eq. 
(6) below. 

Let us describe the essentials (for our purpose) of the B-T model, which is a 
way to implement exact Poincare group transformations for a finite, fixed, number of 
interacting particles. In order to define the Poincare transformations, one introduces 
another set of variables, namely, total momentum P, internal momenta ki,...k n 
(2 ki = 0) and internal spins S^. The unitary transformation which relates the 
previous wave functions to the wave functions _ s (P, k 2 , ■ ■ ■ k n ), depending on 
the internal variables, is the following: 



* sl ,... Sn (pi,...p„) 



\ M 



:n^4) e (3) 
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(]jD i (R i ) SitS 0¥^ s ,jZp i ,k 2 ,...k r , 



i=i 



where (on the right), the vectors ki, the O-components k® and p®, M , and the Wigner 
rotations Ri are functions of the Pi defined as follows : 

P" = + ml M = ^(Ep,-) 2 , (4) 

(notations : B p is the boost (Vp 2 ", 0) — * p, -Dj(-R) is the matrix of the rotation R for 
the spin Si). Let us stress some virtues of unitarity: starting from an orthonormal 
set of internal wave functions, one gets an orthonormal set of wave functions in any 
frame. 

The Poincare generators are then defined as P for the space translations and, 



H = P° = \/P 2 + M 2 (5) 
(9 n n 8 

j = -iPx^ + s, s = J2s' i -i^2k i x 

OP i= i i=2 dh 



where, in order to satisfy the Poincare commutators, the sole requirement is that the 
mass operator M depends only on the internal variables and is invariant by rotation. 
Namely, M must commute with P, ^jg and S. It is important to realize that, in the 
interacting case, we have to deal simultaneously with two different mass operators, 
M of equation (4), and the true mass operator M which appears in (5) and contains 
the interaction. 



The non-interacting case corresponds to M — M = X)"=i V ^1 + m i and, 
through the transformation (3), the generators (5) reduce then to the sum of the 
one-particle free generators, which is the main virtue of (3). 

It is useful to notice that, if M stands for the mass of a particle and S for its 
spin, formulae (5) gives precisely the one-particle free generators when (as assumed 
here) the moving spin states are defined by 

\P,fi) = J^B P \6,ti) (6) 

(P° = ^P 2 + M 2 , (P\ rf\P, p) = (2tt) 3 6(P> - P) 8^) 
The generators follows from the finite transformations, which are given by : 



A|P,/x> = J { -^J2 d ( b ap ab p)»>,» |A*V> (7) 
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where D(R) is the matrix of the rotation R for the spin S. 

First it is then easy to see why (5) satisfy the Poincare Lie algebra, since the 
calculation of the commutators in (5) is the same as in the one-particle free case, due 
to the commutativity of the mass operator M with P, -Q* and S. 

Next the finite transformations generated by (5) are just given (on eigenstates 
of M and P) by (7), except that, instead of acting on the spin S, the Wigner rotation 
B^pABp applies now to the internal spins S[ and the internal momenta k; L . In fact 
they are directly given by (7) on eigenstates of M, P, S 2 , S z . 

We are now in position to construct the wave functions of moving bound states. 
Let v ? si,...s„(^2, • • • k n ) be an eigenstate of M, S 2 , S z (in the Hilbert space reduced with 
respect to P). The corresponding (generalized) eigenstate of M, P, S 2 , S z with P = 
writes : 

<:2f'(£, . . . k n ) = (2tt) 3 S(Q) <p sl ,... Sn (h, ...k n ) (8) 
The moving wave function v|>( p )' m * is obtained by applying the boost Bp (P° = 



P 2 + M 2 ) to the wave function (8). Since we start from P = 0, we may use eq. 
(6), and the result is simply : 

k 2 , . . . k n ) = (27T) 3 5(Q ~ P) <P sl ,...s n (k2, ...K) (9) 

The wave function in the one-particle variables, for a bound state of momentum P, 
is then obtained from (3) : 



Erf » Jkl 



\w M^ < 10 > 



J_ 

1=1 



*£.. a „(pi>---&) 

n 

E (II AW^) (2n) 3 5m-P)ip s[ ,... s , n (k 2 ,...k n ) 



When P — 0, this reduces to : 

*2... S „(Pl, • • -Pn) = (2n) 3 6(m Vs u ..,M, ■ --Pn) (11) 

showing that ip is just the rest-frame internal bound state wave function. 

Introducing the momentum eigenstates given by (10) in formula (1), one gets: 



P \°\ p ) - 1 (11 (ar) ,) ^ MqM % ^ (12) 



E E ■■■*£) 

si, "- s " 

[^(^-^OCp^pOAC^i)],;^ [ft A(ir^kJ V.i,..^.(*2, ••• 4) 



i=2 
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where (on the right) the quantities fcj, M°, -Rj are the functions of p 2 . . .p n given 
by (4) with p\ = P — J2i= 2 Pi > an d the analogous primed quantities are also given by 

— * 

(4) replacing p\ by p[ — P' — Y^l= 2 Pi and m\ by m' 1 . 

As said before, eq. (12) does not give a in general a covariant model for the 
current matrix elements. The reason is that the one-particle operator O is covariant 
with respect to the free one-particle Lorentz transformations, while the transforma- 
tions (5) depend on the interaction. However, it turns out that (12) becomes covariant 
in the heavy mass limit mi,m[ — > oo. This is a general result which requires only the 
covariance of O for free one-particle transformations. We leave the proof to a later 
publication and restrict in this letter to the case of mesons, namely systems of two 
spin 1/2 particles. Furthermore, in accordance with heavy mass limit ideas in QCD, 
we assume spin-independent interaction, and satisfy in this case Isgur-Wise scaling. 
We consider only the ground states, pseudoscalar and vector mesons. 



3 Formulation in terms of Dirac matrices. 

Under the assumption of spin independent forces, the pseudoscalar and vector wave 
functions write 

<Ps lt 82(k 2 ) = -j= (0- 2 ) S1 ,S2 <P(k 2 ) (13) 

where ip(k 2 ) is only required to be invariant by rotation, and e is the rest frame 
polarisation vector of the vector meson. The spin sums in eq. (12) reduce to the 
trace of a 2 x 2 matrix. For example, an harmonic oscillator potential leads to: 

v {k 2 ) = (2tt)§ (R/V^f e-#%' 2 
Using the relation a 2 D(R)a2 = D(R^ 1 ) t , one gets: 



(P'e'\0\Pe)= f-™- ' ir l] ^ ' ! ^ 2 2 (U) 
{ ' ' ' J (2tt) 3 \ M>M ^ p\ U4j 



fpjvj 

l - Tr[(e'. ( ?)tD( J Rr 1 )0(p 1 ',p 1 )D( J R 1 )(e. ( x)D( J R 2 1 J R' 2 )] y,'(k 2 )*y?(k 2 ) 

for the matrix element between vector mesons. The other matrix elements are ob- 
tained by omitting (e.a) or (e'.a) or both under the trace. 

Next, this formula can be written in a more familiar form, involving the 4x4 
Dirac matrices instead of the 2x2 Pauli matrices. Between vector mesons, one gets 



8 



16 



Tr[0( mi + + thM™2 + fc)7fc#(l + f)« + #)] ^'(k 2 ) V(k 2 ) 



u'°u° 



Pi Pi ijk? + mi ^ + mi ^ + m 2 + m 2 
In (15), the following additional notations are used. The unit 4- vectors u and u' : 

Pi + Vi 



(16) 



P1+P2 



17) 



The 4-vectors e M and e^/ are related to the polarisation 4-vectors e = Bp(0,e) and 
e' = B P ,(0,e") by 

£u = B u B P e, e u , = B u iB P ,t (18) 

Moreover, the O under the trace stands for the Dirac matrix appropriate to the 
current considered, for example 7^, 7^75, etc. The matrix elements for other mesons 
are obtained by omitting 7 5 / u or 75/^* or both under the trace. 

Let us describe the main steps to deduce (15, 16) from (14). 0(p',p) s ' tS is the 
matrix element of 7°0 between spinors of the form 



— R ( Xs 

p° p V 



lp° + m 



/ Xs \ 



2p° 



p. a 



V p°+m^ s ) 

where the boost B v also stands for its matrix in the Dirac representation : 

m + fry 

B„ = 



(19) 



^2m(p° + m) 



(20) 



The 2x2 Pauli matrix in (14) may be considered as a 4 x 4 Dirac matrix with only 
its upper left 2x2 block non vanishing. The matrix 0(p',p) then writes : 



/mjmi l+7° p -i^ D 1+7 C 
0(Pi,Pi) = ; — z — B„/ OB, 



Vi ^^pi 2 



(21) 



(using 7°B P 7° = B~ l ). Also we replace the Wigner rotations in (14) by their ex- 
pressions (4) as products of three boosts. The resulting expressions contain matrices 
sandwiched between B u and B~ x and between B u / and B~}. They are reduced using 



B„(l + /)B i -'B„- 1 = (l + ?i) . m ^/ 2 ; 

^2m 2 {k% + m 2y 



(22) 
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and similar other formulae. (22) is obtained from the relation BufcB^ 1 = ^.BuX, 
which simply expresses that the 7^ matrices are forming a 4-vector, since from (20) 
we have 

(1+7°)^ = (1 + 7°) , m ) + h = (23) 

y2m 2 (fc 2 + m 2 ) 

and u = B u (l, 0), p 2 = B u k 2 . Finally, we use 

B^e.a^B- 1 = 75^ B u> (e'*a)^ B~} = 75 # (24) 



4 Heavy mass limit: covariance and Isgur-Wise 
scaling. 



We now consider the heavy mass limit of (15, 16). This limit is defined as m 1? m[ — > 00 
£r and v = £ 



with v' — j^r and v — ■£ fixed (M' and M are the final and initial meson masses). 



It is also assumed that — — > 1, ^ — > 1. It is then found that the integrand in (15) 



has a limit for fixed integration variable p 2 . We have 

Pi P'1 



, V, (25) 

mi m x 

it — > i>, 1/ — > 1/ 

1,0 jL/0 
fc l -1 I 
> / 

mi m x 

A; 2 - S: 1 ^, ^ - B~}p 2 



Also, 



(B-V2) = pa.t;, WW = P2V (26) 

by invariance of the scalar product since B~ l v = (1,0), B~}v' = (1,0). The limit of 
(15, 16) is therefore 

(P e'\0\P e) = [ ^ 2 1 (2?) 

2vV^i (^PZyJfa.v' + mJ&.v + mJ 

1 



- Tr[0 75 /(1 + f)d>2 + m 2 )(l + y)75<*] ^(B^f^B^p,) 

Now, this expression is explicitly a covariant function of the 4-vectors P' and P (or 
i> ' and v) because we have the integral with the invariant measure ^ of a covariant 
function of p 2 , v' and i>. The only point perhaps not immediately apparent is the in- 
variance of Lp'(B~i l p 2 )*Lp(B~ l p 2 ). In fact, ip'(k) and (p(k) being invariant by rotation 
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are function of k 2 and, according to (26), we have 



{B-h2? = (P2-v') 2 - ml, (B~ 1 p 2 ) 2 = ( P2 .v) 2 - ml (28) 



Therefore p'(B~, 1 p2)*(p(B~ 1 p2) is a Lorentz scalar. 

Using its covariance properties, eq. (27) can be reduced further. Indeed, (27) 
can apparently be expressed (for all Dirac matrices O) in term of three independent 
form factors A, B, B', defined by : 



A (v -v) = I ^7= =, w ; = V ( B v P2) <P(B V p 2 ) (29) 
P2 J{p 2 .v' + m 2 ){p2-v + m 2 ) 



B\v'.v)v^ + B(v / .v)v" 

J l Z7r J P2 ^(p 2 .v' + m 2 ){P2-v + m 2 ) 

After integration however, the expression (1 + f)(j> 2 + m 2 )(l + becomes 

(1 + f){B'f + B^ + m 2 A){\ + f) = (B + B' + m 2 A)(l + ^)(1 + ^') 

and, in fact, only the combination £ = + + m 2 A appears. We obtain finally the 
standard scaling formula : 

(P',e'\0\P,e) = \ Tr[0 75 /(1 + f)(l + /Wl £(v'.v) (30) 

with the Isgur-Wise function £ given by : 



>/ , x J f dp 2 y/(p2.v')(p2.v) 

p 2 .(v' + v) + m 2 (v'.v + 1) , 



^ {p 2 .v' + m 2 )(p 2 .v + m 2 ) 



V \B-}p2f^{B- 1 p 2 ) 



The matrix element with one or two pseudoscalar mesons are obtained from (30) by 
omitting 75/ or 75/'* or both under the trace. 
One may notice that 

iiX) = / Mp 2 ) (32) 

and that flavor independent forces and heavy mass limit entail <// = ip, so that 

an 1. 

Although they do not give a general expression like (27) and provide their 
results only in a particular frame, by gathering the various indications given in the 
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paper by Close and Wambach [p!5[ , we apparently end up with precisely this expres- 
sion; we have also checked that their expression of the slope p 2 , given for harmonic 
oscillator wave functions, coincides with ours. An important advantage of our ap- 
proach is that we derive it from the general Bakamjian-Thomas formalism and we 
demonstrate the covariance and scaling properties of the result. 

Let us write the vector [O = V 1 = 7^) and axial (O = A^ 1 = 7^75) matrix 
elements with an initial pseudoscalar, following from (30) : 

(p'\V^\P) = - (v ,fl + v") £(v'.v) (33) 



vW (P',e'\V»\P) = -- 5>^>p4* av'-v) 



vW (P'\A^\P) = 

vVM (P', e'\A»\P) = l - [{y'.v + 1) - (v.e'*) v'»] £(v'.v) 



The corresponding form factors are 



f + (q 2 ) = V(q 2 ) = A 2 (q 2 ) = A (q 2 ) = ^=t£ (34) 

^^^^^^^ 

(v .v = ) 

v 2M'M ' 

5 Lower bound on the p 2 slope. 

Let us now consider the slope p 2 = — £'(1) of the Isgur-Wise function given by (31). 
We establish here the following optimal lower bound of p 2 : 

P 2 > \ (35) 

Eq. (31) is of the form 



with 



£(v'.v) = / — F(p.v',p.v,v'.v) (36) 
J p 



7-1/ / \ 1 \J x'x x' + X + m2(?/ + 1) , , 2 2\* j?/ 2 2\ 

F{x,x,y) = —— f(x -m 2 ) f(x -m 2 ) 

fik 2 ) = V (k) (37) 
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The integral in (36) depends only on the scalar product v'.v because the integration 
measure dp/p° is Lorentz invariant. Expanding (36) around v'.v = 1, we find 



£'(1) = -]% [5 J? 2 WVV, 1) " &*W, 1)] (38) 



dp r l ^ 2 

And using eq. (37) for F, the following formula for the slope is obtained 



p2 = llj& [VpV(p)]*-[VpV(p)1 (39) 



f dp 2 1 mj 1 m 2 * 
+ i(W [ 3 + 4W~3^; ] ^ ^ 



(2tt) 3 l 3 4 (p ) 2 3 p°+m 2 1 

It is obvious that p 2 > 0. In fact this expression is a positive quadratic form of the 
wave function ip, and the best lower bound of p 2 is given by the greatest lower bound 
B of the spectrum of the following corresponding self-adjoint operator T: 

r = -3" A ^ + 3 + 4#-3^ (40) 

It is easily seen that the lower bound is obtained in the S-wave subspace. Indeed, 
the reduction T L of T to the L-orbital eigenspace, acting on the radial functions f{p) 
related to (p by ip(p) = \p | _1 /(|p \)Yif{p), writes 

1 »l_ L(L+l)W + 2 1 mi__l^_ 
L 3 P dp 2P 3 p 2 3 4 (p ) 2 3 p°+m 2 1 ; 

and, due to the "centrifugal barrier", we have T L < T L > if L < L' . 

We may therefore concentrate on (41) with L = 0. Numerical integration of the 
ordinary differential equation T f = B'f (with initial condition /(0) = 0, /'(0) = 1) 
gives a strong indication that the lower bound of the spectrum is 

B = \. (42) 

Indeed, one finds that every B' > | is a generalized eigenvalue and that the (non 
normalisable) eigenfunction is oscillating for B' > | and non oscillating for 5' = |. 

For a proof of (42), we use a unitary transformation U : L 2 ([0, 00 [) — > L 2 ([0, 00 [) 
which has the virtue of converting T L into an ordinary Schrodinger operator. U is 
just the following change of variable : 



(Uf)(x) = A/m 2 cosh(x) /(m 2 sinh(x)), (43) 



(U- 1 f)(p) = ^=f (Arccosh^- ) 



p° " m 2 
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The operator p°d/dp becomes : 



o ^ TT- 1 ^ ^ 



=^-2 tanh W < 44) 



and the operator T (eq. (41) with L = 0) becomes 



3 1 d 2 1 

[/T ET 1 = ~ + i (--^ - L -—-) (45) 

4 3 v dx 2 1 + cosh(x) ' v 7 



Here we can see that B < | because the (always positive) expectation value of the 
operator —d 2 /dx 2 on a spreading function f\(x) = vXf(Xx), X — > 0, (||,/a|| = 



1) goes to 0. On the other hand we have B > | because the operator between 
parenthesis in (45) is positive, as can be seen from the following identity: 

d 2 1 ,d 1 .., d 



dx 2 1 + cosh(x) y dx sinh(x) ^ dx sinh(x) ^ 

Finally, one sees as follows that the lower bound B = 3/4 of p 2 cannot be 
attained. If p 2 = 3/4, the expectation value of (46) vanishes on the normalized 
function f(x) which corresponds to <p(p). This implies ( — ^ + sin l^ )f(x) = 0. 
However, the solutions f(x) = c tanh(|) of this equation are not normalizable. 



The bound (35) is obviously in agreement with Bjorken's lower bound JT7] 



p 2 > 1/4. It does not contradict Voloshin's upper bound |L8| p 2 < 0.75 ± 0.15, but 



not much room is left and a careful study of the relativistic generalisation of Thomas- 
Reiche-Kuhn sum rule in our model should be performed. Neither does it contradict 
de Rafael-Taron's |T9[ conservative estimate p 2 < 1.7 not to speak of their rigorous 



bound p 2 < 6.0. 

Let us emphasize that the bound 5 = 3/4 applies to a large class of quark 
models, which at least in this heavy quark limit show very remarkable properties, 
and is independent of particular parameters of the models. Whence the interest of 
testing it experimentally. But this is not an easy task. Let us recall that p 2 is not 
directly comparable to the p 2 measured by CLEO II experiment [^TJ, from which 
it should differ by QCD radiative corrections and l/m c corrections PDf. The QCD 



radiative corrections enhance p 2 versus p 2 ; the l/m c corrections are not known from 
exact QCD, and cannot be safely deduced from our model in which they are not 
covariant. The central value of the present CLEO II data: p 2 = 0.87 ± 0.12 ± 0.08 
PH is compatible with our lower bound, although not far from it. But one must keep 
in mind the above-mentioned corrections. A further improvement of data may prove 
very instructive. 

The wave functions that saturate the lower bound are, for example, of the 

type: 

<p{p) oc V~e (p )-( L5 + £ ) (47) 
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for which p 2 — > 0.75 when e — > 0. 

Harmonic oscillator wave functions, 



<p(£) = (2tt)I (i?/v^ 



I e -flV/2 



(48) 



which are most commonly used, yield p 2 well above the bound 0.75. Indeed 
we found that the lower bound for harmonic oscillator wave functions is 1.208. In the 
limit m 2 R 2 = they give p 2 = 5/4. The minimum 1.208 is obtained for m 2 R 2 ~ 0.05 
which is far below the standard quark model parameters, say, m 2 ~ 0.1 GeV 2 , R 2 ~ 6 
GeV -2 , leading to m 2 R 2 ~ 0.6. Using the latter values we obtain p 2 ~ 1.37 which 
is, as expected, very similar to Close and Wambach's [15|: p 2 ~ (1.19) 2 . In the weak 
coupling limit m 2 R 2 — > oo one gets 



m 2 R 2 



+ 1 + 



m 



'R 2 



(49) 



and the Isgur-Wise function takes the very simple form in the large m 2 R 2 : 



av.v') = 2 - L^ eM-™ 2 R 2 (v.v> - l)/2] (l + O (-^s)) (50) 

We have checked numerically that expression (50) is a surprisingly good ap- 
proximation to the Isgur-Wise function obtained with the standard quark model 
parameters: m 2 ~ 0.1 GeV 2 , R 2 ~ 6 GeV -2 , although m 2 R 2 does not look so large, 
this amusingly simple formula is not phenomenologically very useful since, as already 
stressed, the harmonic oscillator wave functions yield a p 2 much above the lower 
bound and above experiment. 
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